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ABSTRACT 

Prior to the explosive burning of a white dwarf (WD) that makes a Type la supernova (SN la), the star 
"simmers" for ~ 10 3 yrs in a convecting, carbon burning region. I estimate the excitation of g-modes by 
convection during this phase and explore their possible affect on the WD. As these modes propagate from 
the core of the WD toward its surface, their amplitudes grow with decreasing density. Once the modes reach 
nonlinear amplitudes, they break and deposit their energy into a shell of mass ~ 10~ 4 M Q . This raises the surface 
temperature by m 6 x 10 8 K, which is sufficient to ignite a layer of helium, as is expected to exist for some 
SN la scenarios. This predominantly synthesizes 28 Si, 32 S, 40 Ca, and some 44 Ti. These ashes are expanded 
out with the subsequent explosion up to velocities of ~ 20,000 km s" 1 , which may explain the high velocity 
features (HVFs) seen in many SNe la. The appearance of HVFs would therefore be a useful discriminant for 
determining between progenitors, since a flammable helium-rich layer will not be present for accretion from a 
C/O WD as in a merger scenario. I also discuss the implications of 44 Ti production. 
Subject headings: convection — stars: oscillations — supemovae: general — white dwarfs 



1. INTRODUCTION 

The use of Type la supernovae (SNe la) as cosmological 
distance indicators has brought attention to the uncertainties 
that remain about these events. It is generally agreed that they 
result from the unstable thermonuclear ignition of a C/O white 
dwarf (WD), but the exact progenitor is still unclear. It may 
be a si ngle degenerate (WD a ccreting from an evolved com- 
panion; Whelan & Ib enll973[). double degenerate (merging of 
twoWDs nWebbinklll984t llben & Tulukovl [19841: iPaczvnskil 
1985), or combination of both scenarios. Any observational 
or theoretical clues that could help unravel this mystery are 
extremely useful. 

In cases when the WD first unstably ignites carbon at 
its center, it subsequently undergoes ~ 10 3 yrs of convec- 
tive simmering before the explosive burning w ave is born 
dWooslev et al.ll2004t IWunsch & Wooslevl 120041) . In single 
degenerate scenarios, central ignition occurs when the accre- 
tion proceeds at rates slower than the WD thermal timescale 
dHernanz et al.| [l988). Central ignition can also occur via fo- 
cusing of a shock from a surface helium det onation, but thi s 
leads to core detonation and not convection dFink et al.l2010t) . 
During the late stages of convection there is considerable 
luminosity going into convective motions (~ 10 45 ergs s -1 ), 
which is much greater than the Eddington luminosity of a 
Chandrasekhar WD of Lem ~ 2 x 10 38 ergs s" 1 . This energy 
is expected to be bottled up within the convective region be- 
cause the thermal conduction timescale is ~ 10 6 yrs, which 
is much longer than the convective timescale. If just a small 
fraction of this energy could be transported closer to the WD 
surface, it might have an effect on the surface structure, and 
may even have important observable consequences. 

In the present work I consider the stochastic excitation of 
g-modes by convection and how they may affect the WD sur- 
face. The presence of such modes was first suggested by 
iPiro & Chand (12008ft and subse quently observed in the sim- 
ulations of Zingale et al. (2009). In ^] I make analytic es- 
timates for the luminosity and total integrated energy in g- 
modes. I show how the modes grow as they propagate to- 



ward shallower densities, and argue that the mode energy is 
deposited near the WD surface when the modes break due 
to reaching nonlinear amplitudes. In §|3] I study the detailed 
structure of the heated WD surface layers. The temperature 
rises sufficiently to ignite a surface helium layer, and the re- 
sulting ashes are predominantly composed of 28 Si, 32 S, 40 Ca, 
and per haps 44 Ti. This may e xplain the high velocity features 
(HVFs, iMazzali et al]|2005bl) seen in many SNe la. In gj I 
summarize this study and discuss where possible future work 
is needed. 

2. G-MODE EXCITATION AND PROPAGATION 

Ignition of 12 C occurs when th e heating from carbo n fu- 
sion overpowers neutrino cooling (Nom oto et al.lll984b . The 
central temperature T c then rises and a convective core grows, 
eventually encompassing ~ 1M Q of the WD after ~ 10 3 yrs. 
For a central de nsity p c , the energy generation rate from car- 
bon burning is ( Wooslev et al. 20QJ) 
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where r c .g = r c /10 8 K and p c g = p c /l0 9 g cm" 3 , and equal 
mass fractions of carbon and oxygen are assumed. The 
central temperature increases on a heating timescale ?/, = 
(d\nT c I 'dt)~ l , which gets shorter as T c becomes larger, and 
generally depends o n the size of the convecting region 
(IPiro & Changl 120081) . The simmering ends once f/, < t c , 
where t c is the eddy overturn timescale. At these late times, 
individual e ddies may experience significa nt heating during 
their transit dGarcia-Senz & Wooslevlll995l) . and there is not 
sufficient time for the entire convective region to respond to 
the increasing T c . In this case it can be approximated that 
th ~ CpT c /e, where c p is the specific heat capacity at con- 
stant pressure. The heat capacity in the WD core, including 
Coulomb corrections, is c p « 1 .3 x 10 7 ergs g -1 K -1 , which is 
used to find 
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Woos ley et al. (2004), using the KEPLER code dWeaver et all 
119781) . estimate that convection ends when T c w 7.8 x 10 8 K 
and p c w 2.6 x 10 9 g cm" 3 . This gives f/, ~ 10 s, roughly the 
convective overturn timescale as shown below. 

Integrating over the burning regi on of the core, the total 
luminosity carried by convection is (Woosle y et al.H2004b 
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The convective velocity at the largest scales is V c « 
(L/4Trr 2 p) 1 ' 3 . What is crucial for driving g-modes is the prop- 
erties near the top of the convective zone, which has density 
and radius p, and r, , respectively. The velocity here is 
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where p, % = p,/10 8 g cm" 3 and r, % = r,/10 8 cm. The spectrum 
of g-modes excited by the convection is peaked at a frequency 
equal to the eddy turnover frequency lu c w V c /H t , where H, 
is the scaleheight at the top of the convection. Taking H, » 

2 x 10 1 p\ g^io cm > where g w = g/\0 10 cm s" 2 , results in 
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Note the convective timescale is t c ~ H,/V c ~ 10 s, roughly in 
agreement with when the convection should end, as discussed 
above. These waves propagate in the non-convective WD sur- 
face layers if their frequency satisfies uj c < N, where /V is the 
Brunt-Vaisala frequency. This is approximated as 
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where p (no subscript) is the density at some position near the 
WD surface, p% = p/10 8 g cm" 3 , is Boltzmann's constant, 
T, is the temperature at the top of the convection, Z is the aver- 
age charge per ion, and Ep is the Fermi energy for a degener- 
ate, relativistic electron gas. I take Z = 13.8, as is appropriate 
for a mixture of equal parts carbon and oxygen, and ignore the 
scalings with composition to simplify the presentation. Since 
uj c < N at the convective boundary, and at shallower depths 
N oc /0~ 1//3 , the g-modes propagate freely toward the surface. 

The fraction of L c that can be put into g-modes is di- 
rectly proportional to the Mach number of the convective ed- 
dies near the top of the convective zone (Goldrei ch & Kumarl 
1990fl For a soundspeed c s = (4P/3p)^ 2 , 
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The g-mode luminosity is then 
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' Some energy is expected to go into />modes as well, but because the 
efficiency is oc Ma' 5 / 2 , this is negligible. 



Comparing the dependence on T c in equations © and (O 
shows that L g oc f^ 1 4 . The total amount of energy put into 
g-modes up to any given time is therefore 
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Depending the final T c , about ~ 10 46 ergs goes into g-modes. 

As the g-modes propagate into the non-convective surface 
layers, they satisfy the dispersion relation 
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and have a group velocity of V g = uj c /k r , where k r and kh are 
the radial and horizontal wavenumbers, respectively. This re- 
lation does not include rotational modifications. The radial 
wavenumber is unaffected by the Corio lis force if the spin £7 
is small in comparison to the buoyancy ( Chapman & Lindzen 
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In this limit the mode equations can be simplified using the 
"traditional approximation" to separate vertical and horizon- 
tal parts. Although k r remains the same, kh can depend on the 
angle with respect to the rotation axis, since the Coriolis force 
pushes the modes to be more concentrated near the equa- 
tor (s ee the angular eigenfunctions plotted in iPiro & Bildstenl 
2004). The spin of accreting WDs may indeed exceed that 
given by equation ( fTTT i, but rotational modifications are ig- 
nored to simplify the current study. 

The total luminosity in g-modes, L g , is carried by n modes, 
each with a characteristic Lagrangian displacement £. As- 
suming that the modes are excited with random phase, the 
luminosity is related to the size of the perturbations via 



L g « 4irr 2 p(ujO 2 nV g 
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From incompressibility, the components of £ = £ 2 + £ ; 2 are 
related by (^ r /^/,) 2 « (kh/k r ) 2 . The linearity of the modes 
is best represented by the dimensionless quantity k r ^ r , since 
both density inversion instabilities and Kelvin-Helmholtz in - 
stability set in when k r ^ r — 0(1) dMihalas & ToomrelfT98Tl) . 
Using the above relations, 
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The number of modes is found by integrating over the mode 
spectrum and the surface area of the convection. There re- 
main uncertainties about t he power spectrum of the convec- 
tion ( see the discussions in lWooslev et a l. 2004; Kuhlen et al.l 
2006), which in turn determines the spectrum of modes. It is 
therefore sufficient to estimate 
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keeping in mind that this is a lower limit since a wider spec- 
trum will require more modes to carry the same luminosity. 
In the limit uj <C N, the amplitude simplifies to 
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and taking kf, ~ 1 /H t , I evaluate this as a function of p, 

Ki r « 3 x iQ- 2 g-ilp)T 2 nXi A pf n - (16) 

As the mode propagates toward the surface oc p~ x / 2 N l l 2 oc 
p" 2 / 3 , so this dimensionless amplitude grows. This quantity is 
independent of the strength of convection, as apparent from 
the lack of a dependence on T c or p c . This is because the 

amplitude is k, ^ r ccL g ' u> . As the driving becomes more 
vigorous, the g-mode luminosity increases, but so does the 
frequency, and the two effects exactly balancfl 

Near the surface the thermal conduction timescale becomes 
shorter. The thermal conductivity from electron-ion scattering 
is (Yakovlev & Urpin 1980) 



where m* w Ep/c 2 , n e is the number density of electrons, 
and v = 4m*Ze 4 A / 4-kJi 3 is the electron-ion collision frequency 
with A w 1 is a weakly dependent function of density in the 
ocean and h is Planck's constant. For a lengthscale A, the local 
thermal time is t± ~ pc p X 2 /K, giving 

f ^ 7xl ° 4 ^ /3 (^ik^) 2yrs ' (18) 

where A has been scaled to the typical horizontal wavelength 
of the g-modes. This is much longer than the timescale for 
a mode to travel a scaleheight t g w H/V g « (H / H t ~)N / u) c ~ 
40 s, so conduction does not damp the modes. Since breaking 
occurs before damping, setting k r ^ r w 1, 1 find 

c , n 5 -3/2 11/8 -3/2™3/8 -3 , 1ri s 

for the mode breaking density. 

3. SURFACE HEATING 

For the breaking depth found in equation (fT9] l, the corre- 
sponding mass of material is Mb ~ 4irr 2 pbH(pb) ~ 10~ 4 M . 
If an energy E g is put into a shell with mass Mb, it does 
not eject the material, since the binding energy for a Chan- 
drasekhar WD is much greater at GMM\,jR ~ 10 47 ergs. In- 
stead the energy input heats the surface, which I now explore. 

A simple estimate for the change of temperature is AT ~ 
E g /cpMb ~ 10 9 K, but this does not take into account how the 
heating is distributed. The heating at a depth pb acts like a hot 
plate, which still cannot be carried by conduction due to the 
long thermal timescale (see eq. [1181 . evaluated at p/,). Instead, 
a secondary convective zone grows, and since the heating is 
now distributed over this entire new convective region, the 
temperature rise at the base will in general be less than the 
AT ~ 10 9 K estimated above. 

To better understand the extent and peak temperature of this 
secondary convection zone, I build a series hydrostatic surface 
models where the convection zone is treated as an adiabat. 
The initial model is assumed isothermal with a temperature of 
Ti = 6 x 10 7 K, similar to the surfac e temperature in accreting 
WD models dYoon & Langer 2004). The convection extends 
from a bottom mass coordinate M\ set by the g-mode breaking 

2 In detail, as the convection proceeds, r t will increase, while g and p, 
will decrease, but these are small correction with respect to the simplifying 
assumptions I am making here. 
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FIG. 1. — Example temperature profiles for the surface convection zone 
caused by energy deposited from g-modes. In the top panel I use M = 
\31Mq,R = 1.6 X 10 s cm, and pi, = 3 X 10 5 g cm 4 , and in the bottom panel 
M = 1.35M , R = 2.1 X 10 8 cm, and p b = 10 5 g cnT 3 . The thick lines show 
convective profiles for different amounts of total internal energy, as labeled. 
The thin, short-dashed lines are curves of constant ty a , demonstrating when 
helium burning begins. 

depth up to top mass coordinate Mi set by the total energy 
injected by the g-modes 

E g = / c p [T(M r )-Ti]dM n (20) 

where T(M r ) is the temperature profile. This integral is fairly 
insensitive to the assumption of 7), since at late times the con- 
vection has much more energy than the initial temperature 
profile. Even though Mi is set fixed by the mode breaking 
depth, the density at the base of the convective zone decreases 
because of the increasing temperature at fixed pressure. 

In Figure Q] I plot example thermal profiles that would be 
present from energy deposited by g-modes. In the top panel I 
use M = 1 31 M Q , R = 1 .6 x 10 s cm, and p h = 3 x 10 5 g cm -3 , 
and in the bottom panel M = 1.35M , R = 2.1 x 10 8 cm, 
and pt, = 10 5 g cm" 3 . The central densities of the two mod- 
els are 6.1 x 10 9 g cm" 3 and 2.0 x 10 9 g cm" 3 , respectively. 
The difference in breaking depths wa s estimated from core- 
convective models (as calculated in iPiro & Bildstenl 120081; 
iPirol 12008). Each thick line shows the thermal profile for 
a labeled amount of g-mode energy E g , solved using equa- 
tion (f20b . The more massive WD is plotted with a larger E g 
due to its higher central density (eq. [|9])- The base tem- 
perature rises by AT ss 6 x 10 8 K and for the top panel and 
w 3 x 10 8 K for the bottom panel, but for a given E g , AT is 
roughly the same. These detailed profiles also provide an es- 
timate of the position of the top of the convection zone, which 
is near ~ few x 10 3 g cm" 3 . If the thermal time here is suffi- 
ciently short, heat leaves the top of the conve ction and stunts 
its gro wth, as is found for Type I X-ray bursts (Weinberg et al. 
2006). Since the thermal time t± ~ 1 yr (eq. 1181 ) is much 
longer than f/,, conduction does not truncate the secondary 
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convective zone. 

Besides changing the entropy and thermal profile, the heat- 
ing may have an additional effect if the surface layers are 
rich in helium. A thin layer of helium is expected in many 
single degenerate scenarios that require accretion to grow to 
a Chandrasekhar mass. For an energy generation rate from 
triple-a reactions of £3 Q , the timescale for increasing the tem- 
perature is ?3 a ttCpT/c3 a . Figure Q] plots ?3 Q = 10, 10 2 , and 
10 3 s (thin, short-dashe d lines), using the reaction rates from 
Fushiki & Lamb1( ll987l) . Comparison of these curves with the 
convective profiles demonstrates that triple-a reactions can 
begin burning the layer. Helium burning favors deeper mode 
breaking and more massive WDs, since a higher density leads 
to a shorter t3 a , A more detailed study is required to assess 
the full range of WD masses and breaking depths needed for 
helium ignition, which is outside the scope of this work. 

The low s hell mass is below wh at is needed for dynam- 
ical burning (IShen & Bild sten 2009]), so a detonation is not 
expected. Nevertheless, energy from helium-burning pro- 
vides fx 5.8 x 10 17 ergs g" 1 , similar to the binding energy 
of a Chandrasekhar WD (~ 10 18 ergsg" 1 ), so it is possi- 
ble that some material is expelled. A breaking depth of 
Pb ~ (3-5) x 10 5 g cm" 3 corresponds to a base pressure of 
~ (4-9) x 10 21 ergs cm" 3 . The burning of helium at constant 
pressure was explored by Hashim oto et al.l ([1983), the prod- 
ucts of which are summarized in their Figure 10. For this 
pressure range, the predominant elements are 28 Si, 32 S, and 
40 Ca, but up to a mass fraction of ~ 0.1 of 36 Ar and 44 Ti may 
also be present depending on pi,. 

The composition and position of these elements are sim- 
ilar to what is required for the HVFs seen in many (or 
perhaps most) SNe la dMazzali et al.l l2005bl) . HVFs most 
prominently show 40 Ca, although 28 Si is sometimes present. 
Currently expla nations for these fea tures include circumstel- 
lar interactions dGerardy e t al. 2004 ]), three-dimensional d en- 
sity or composition enhancements dMazzali et al.ll200"5al) . or 
ashes on the surface from a gravitationally confined explo- 
sion dKasen & Plewal 120051) . Modeling of th e early time 
lightcurve, including the HVFs was done by Tana ka et al.1 
(2008). Surface burning from g-mode energy input suggests 
that 32 S is expected to be fairly abundant as well, but the 
physics of its line creation may not be optimal for seeing it. 
Nevertheless, searching for signs of 32 S would be an impor- 
tant test of this hypothesis. If this model is correct, the pres- 
ence of HVFs would be an indication that helium accretion 
(or hydrogen accretion which then burns to form helium) is 
responsible for these SNe la, which is expected for a single 
degenerate progenitor or a degenerate helium donor (like for 
AM CVn stars). 



4. CONCLUSION AND DISCUSSION 

I considered the driving of g-modes via convection dur- 
ing the simmering stage prior to SNe la. I found that the 
g-modes are driven with a luminosity L g > 10 42 ergs s" , in- 
jecting E g ~ 10 46 ergs at a density of pt, ~ few x 10 5 g cm" 3 
where the modes break. I considered the effect of this en- 
ergy in generating a secondary, surface convective region and 
showed it reaches base temperatures of«6x 10 8 K, sufficient 
to ignite helium. The ashes expected from this burning are 
similar to the elements observed in HVFs, which may mean 
that HVFs indicate a progenitor channel different from the 
merger of two C/O WDs since these would not have helium 
present. 

This study motivates further, more detailed investigations. 
The mode excitation estimates employ ed here make use of 
studie s done in the context of the Sun dGoldreich & Kumarl 
1990), and may need to be modified for a degenerate equa- 
tion of state. In addition, rotation can strongly influ ence the 
morphology of WD convection dKuhlen et alJ 120 06). which 
would in turn alter the mode wavelengths and spectrum (also 
see the discussion of rotational modifications to the modes in 
§2). Rotation also breaks the spherical symmetry of the WD 
and would imprint asymmetries into the flux of g-modes. 

A study of thin helium-shell burning on massive WDs 
would provide detailed predictions for the burning products 
and determine what fraction of material can be expelled. The 
production of 44 Ti is especially interesting for testing this 
model and investigating SNe la progenitors. Our fiducial cal- 
culations estimate ~ 10~ 5 M Q of 44 Ti is pr oduced, although 
this varies strongly as a function of pi,. Borkows ki et al.l 
d2Ol0h find (1 -7) x 1O" 5 M of ^Ti in the remnant G 1.9+0.3 
from the decay product 44 Sc. The 44 Sc is most abundant in the 
remnant's north rim, which is surprising if the 44 Ti is expected 
to be synthesized in the neigh borhood of mo re centrally lo- 
cated Fe-peak elements (Iwamoto et al.lll999h . Further stud- 
ies of the distribution of 44 Sc would help determine whether 
surface helium burning is a viable alternative. Detections of 
the 68 keV decay line from 44 Ti in SN e la by future satel- 
lites like NuSTAR (Harri son et al.l l2010) may also provide an 
important constraint. 
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